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Ограничение математического моделирования в геофизике. Адептью функционально-аналитического направ- 
ления геофизики разрабатьтвают методьт и алгоритмьт в рамках известньх интегро-дифференциальньх уравнений 
математической физики, описьвающих поведение геофизических полей. Их усилия направлень на определениєе не- 
которьїх козффициентов зтих уравнений и внесение соответствующих поправок за "сложность" средьт. Зти козффи- 
циентьт определяются по помощи итерационньх схем, уложенньїх в принципьт Тихонова решения обратньтх задач 
математической физики |. 

Но определение зтих козффициентов по функциям, образующим соответствующие системьт уравнений -- принци- 
пиально неоднозначно. Зто явно очевидно в массовьтх машинньтх расчетах на системах большой размерности |2|. А 
при найменьших погрешностях во входньгх данньїх решения еще и существенно неустойчивь. Одолеть зту неустой- 
чивость на порядок сложнее, чем получить формальньте решения задач. Вдобавок, след входньх функций, по кото- 
рьтюм определяют их козффициенть - двумерен, а восстановить нужно трехмерное поле, по значениям которого най- 
ти козффициентьт при третьей координате, не совпадающей с размерностью следа. 

Формальное решение зтой проблемь дает псевдопространственное распределение, "растворяющееся" среди зк- 
вивалентньх решений задачи. К тому же, займствуя аппарат математической физики, геофизики одновременно уна- 
следовали с ним сугубо евклидовую структуру пространств функций, которая для региональньх построений недо- 
статочна |3). Но ограниченность интегро-дифференциальньх операторов, аппроксимирующих математическую мо- 
дель средьт - лишь методическая сторона проблеми. 

Входнвье данньтве - зто, в общем, двумернье проекции геофизических полей на поверхность Земли, измереннье с 
погрешностями на разреженной нерегулярной сети наблюдений. Они априори получень!т с неизвестньми погрешно- 
стями, позтому в дополнение к сказанньїм вьше методическим проблемам усложняют поведение (и определение) 
функций, входящих в соответствующие модельньт?е уравнения. А развитие компьютерного моделирования привело к 
применению дискретной шкальт входньх функций. Учитьвая, что до сих пор не освоєно прямое получениєе вход- 
ньх данньх для обратньх задач геофизики из баз данньїх пунктов наблюдений (4, перевод непрерьтввньтх входньтх 
данньх в дискретную форму в цифровой картографиий неизбежно вносит дополнительньюе искажения, кроме погреш- 
ностей измерений. Зти проблемьі являются следствиєм методологий получения данньїжх прошлого столетия, которая 
отживаєт свое. Мьт их назьтвваєм методологическими. 

Новьг6й геофизико-математический аппарат. Указанньг6"е вьше, и другие математические и методические проб- 
лемь, обусловленнье неадекватностью аппарата математического моделирования геофизических данньтх, заставля- 
ют геофизиков разрабатьгтвать собственньй геофизико-математический аппарат. В его рамках существенно перера- 
ботаньт известньгюе математические методьт с учетом фактов, присущих сугубо геофизическим явлениям. Начал раз- 
работку своеобразного "геофизического диалекта" В.Н. Страхов |3| в потенциальньїх полях, учтя ортогональность 
сигнала и погрешности и дискретность данньх; продолжил В.В. Гольдин в сейсмике, В.В. Аксенов в злектродина- 
мике и др. 

Вне четкой национальной программь развития технологий потенциальньх полей и создания сети мониторинга 
зтих полей нельзя влиять на реформу принципов получения исходньгх данньгх. Возможно, адекватная нуждам гео- 
физиков-полевиков, и интерпретаторов система непрерьвньїх наблюдений для формирования общедоступной базьт 
данньхх гравимагнитньх полей будет в ближайшей перспективе. Но можно предложить некоторьтве соображения по 
новьтм представлениям математических моделей геологической средьт и геофизических полей в рамки апроксимаци- 
онного подхода |б| к созданию подобньх конструкций. 

Новая постановка обратной задачи гравиметрий. Устойчивьте способьт восстановления потенциала силь тя- 
жести крайне важньт и для детальной гравиметрии и для глобального восстановления фигурьт Земли. В последнее 
время ради повьшшения устойчивости и надежности решений в практических приложениях задач теориий потенциала 
применяют новую величину - градиент потенциала сильт тяжести. Теоретически она бьла определена еще в работе 
І7). Современнье приложения градиента потенциала силь тяжести ориентировань на применение в сочетаний с 
классическим представлением аномалий силь тяжести как вертикальной производной потенциала. 

В то же время значительно более вьтсокую устойчивость и более широкий круг применений (хотя и более слож- 
ное двухступенчатое вьічисление) имеет аналитическое продолжение упомянутого градиента для специально полу- 
ченной в |3| новой аналитической (дифференциальной) модели поля силь тяжести. В ней учтена векторная природа 
силь тяжести, градиент ее изменения в зависимости от уровня кривизнь поверхности измерений. Также определеньт 
погрешности, возникающие при замене прямого оператора известного в аналитическом продолженимй уравнения Пу- 
ассона на уравнениєе сильт тяжести: они пропорциональньт геометрическим параметрам области исследований. 

Вкратце основнье зтапь ее вьтвода иллюстрирует следующая цепочка уравнений: 
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В рамках решения однородной части зтого уравнения возникла так назьтваємая нелинейная граничная задача 
Алексидзе для уравнения Лапласа |8), содержащая в граничньтх условиях значения модуля градиента потенциала си- 
льт тяжести. Ее решение полностью теоретически обосновано для гравитирующих обеектов, ограниченньтх поверх- 
ностями Ляпунова. Определеннье трудности ее решения известньтми численньжми методами обнаженьт в сообщений 
91. 


В общем, для определения модуля градиента сильі тяжести я(х х)а (етади( х) на поверхности оС Земли полу- 
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а поверхностная плотность се ) Є є дО, входящая в вьіражение (2), определяєтся из нелинейного уравнения: 
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Методика решения задачи. Для получения решения уравнения сі в |3| полученьт две итерационньге схемі, к 
сожалению, малопригодньве для практических вьічислений большой размерности. Считаєм, зто происходит в силу 
слабой сходимости соответствующего интеграла (в смьтсле главного значения). А его регуляризация через изобра- 
жение отрезком ряда не гарантирует желаємой точности решения. 

Один из возможньх способов решения уравнения (3) ре бьтть таким. кю его искомое решениє 


с), ЄєдС в виде линейной комбинации с(х)- УУауй) . Тогда его квадрат о )- У арацту(х ри (х ). 
Теперь формализуем зо (3) через введенное обозначениє: 
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В таком изображений задача в ИиТтОоГге СВвОДИТСЯ К классической апроксимационной постановке: определить козф- 
фициентьт дії уравнения (4) так, чтобь вьшолнялось условиє минимизации невязки функционала 


із" (о) )- Ла, а) з тіл (5) 


Отньгюне, смотря какую из популярньтгх в теории обратньх задач метрических норм мьт изберем, задача нахожде- 
ния аналитической аппроксимации (4)-(5) решаєтся в том или ином пространстве с помощью классических вариаци- 
онньх методов, в частности, модификаций метода Ньютона. Решение получаєтся в несколько раз бьгстрее (зависит 
от архитектурь компьютера), а преимущества в точности станут яснь только после испьтаний на полевьтх данньхх. 
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Математические постановки конкретньх геофизических задач сводятся в итоге к определению значений козффици- 
ентов уравнений математической физики, входящих в зти постановки. В русле зтого факта решение задачи Алексидзе в 
виде нелинейного интегрального уравнения сведено к линейной комбинации искомьх решений. Такая альтернативная за- 
дача сведена к классической задаче вариационного исчисления. 
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